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Abstract 

In this search, the chaotic behavior of the Rossler system are analysis by using sensitivities 
depend on initial condition. These are implemented by varying the parameter of system. 

The scope of the work is examined by using the software (Matlab) as programming language to get 
sensitivities depend on initial condition (chaos). 
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1-Introduction 

Rossler systems is introduced in the 1970s as prototype equations with the minimum 
ingredients for continuous times chaos. 

Since the Poincar'e-Bendixson theorem precludes the existence of other than steady 
periodic, attractors in autonomous systems defined in one- or two -dimensional manifolds such as 
the line, the circle, the plane, the sphere, or the torus (Hartman, 1964), the minimal dimension for 
chaos is three. On this basis, Otto Rossler came up with a series of prototype systems of ordinary 
differential equations in three-dimensional phase spaces (Rossler 1976a,c, 1977a, 1979a). Systems He 
also proposed four-dimensional for hyper chaos, that is chaos with more than one positive Lyapunov 
exponent (Rossler 1979a,b). 

Rossler was inspired by the geometry of flows in dimension three and, in particular, by the re- 
injection principle, which is based on the feature of relaxation- type systems to often present a Z-shaped 
slow manifold in their phase space. On this manifold, the motion is slow until an edge is reached 
whereupon the trajectory jumps to the other branch of the manifold, allowing not only for periodic 
relaxation oscillations in dimension two, but also for higher types of relaxation behavior as noted by 
Rossler (1979a). In dimension three, the re-injection can induce chaotic behavior if the motion is 
spiraling out on one branch of manifold). In this way, Rossler invented a series of systems, the most 
famous of which is probably (Rossler 1979a). 

In this section we study the chaotic behavior of Rossler system depend on the definition of 

Gulick which is referred to in section two. 



1 



2- Definitions 

Here many fundamental definitions are used . 

• Definition 1 [Periodic attracting] 

Let x be a periodic -n point for a function f then x is attracting period-n point if x is an 
attracting fixed point of f 1 [Gulick,1992 ] 

• Definition 2 [ lyapunov expoent] 

Let J be abounded interval, and f: J— >J continuously differentiate on J. fix x in J ,and let ^x be 
defined by 

Ax=Lim 1/nlf^l (x) ( 1) 

n— »oo 
provided that the limit exist . in that case tat is the lyapunov exponent of fat x [Gulick,1992] 

• Definition 3 [sensitive dependent on initial condition] 

Let J be an interval, and f:J— >J has a sensitive dependent on initial condition if there exist £>0 
such that for any x G J and any neighborhood N of x ,there exist yGN and n>0 such that 
I f(x)-f(y) I >^ [Deveny,1989] 

• Definition 4 [Chaoic] 

A function f is Chaoic if satisfies at least one of the follwing conditions 
(i) f has apositive lyapunov exponent at each point in its domain 
(ii) f has a sensitive dependent on initial condition on its domain [Gulick,1992] 

• Definition 5 [Capacity and Fractal dimension ] 

Let S be subset of R n , where n=l,2 or 3 the capacity dimension of S is given by Dim C S =Lim 

ln(N(G)/ln(l/G) (2) 

G^O 
If the limit exist and is not integer then S is said to be have Fractal dimension [Gulick,1992] 

• Definition 6 [Bifurcation ] 



Consider the differential equation : 
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x=fn(x) ...(3) 

One is especially concerned how the phase portrait of (3) chang as (i varies ,A value (io where there is 
a basic structural change in this phase portrait is called a bifurcation point [Gulick,1992] 



• Definition 7 [Bifurcation diagram ] 

One method of displaying the points at which a parameterized family of function { f M } bifurcates 
and is designed to give information about the behavior of higher interates of arbitrary member of the 
domain of f M for all value of parameter |i [Gulick,1992] 



3-Rossler Model : 

Rossler was able to obtain the simplest nonlinear vector field capable of generating chaotic 
behavior [Rossler, 1976]see however, [Sprott,1994] This attractor is written in the following form : 



x =-x-y 

y = x+ay (4) 

z'=b+z(x-c) 

such that it has a single nonlinear term xz in z * . By fixing a and b in the value a=b=0.2, one has a 
period-doubling route to chaos where a period-2 orbit is created at c=2.6, and being c~4.2 the 
accumulation point of the period doubling cascade, beyond which one has deterministic chaos, 
excepting for the presence of a number of periodic windows. The system has an unstable fixed point 
near the origin whose 2D unstable manifold presumably spans the strange attractor. It appears that the 
strange attractor does not exhibit a remerging tree (or period-doubling reversal) [Stone, 1993], at least 
for not too large values. 



4- Plots Description : 

In Fig.(l) one can see the scatter-plots for the Rossler attractor. 
The left column of plots Fig. (la,lc and le) are the results for the new algorithm, whereas the 
column on the right-hand side, Fig.( lb, Id and If) shows the results for the Wolf algorithm. Both plots 
la and lb have the x-coordinate of the Rossler attractor as abscissa. Analogously, plots lc and Id have 
the y-coordinate of the Rossler attractor as abscissa and plots le and If the z-coordinate. The ordinate 
of all cases is the value of the positive local Lyapunov exponent ^1 (t). 

Fig. ( 2 ) shows the pair wise Renyi spectra corresponding to the plots of Figs. (1.) The dashed 
line is the spectrum for the Wolf algorithm and the full line for the new one. Specifically, parts Rossler 
a, Rossler b and Rossler c denote for the pairs of spectra that correspond to the pairs of point sets 
(la,lb), (lc,ld), (le,lf), respectively [Grond andDiebner,2005]. 
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Fig. (1) Plots of the x-, y-, and z-coordinates of the Rossler attractor against the local Lyapunov 
exponent kl. The left column (a, c, and e) shows the results for the new algorithm, the right column (b, 
d, and f) for the Wolf algorithm. 
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Fig. 2. Plots of the Renyi spectra computed from the point sets of Figs. 1 and 2. The left column shows 
the results for the Rossler attractor, The curve corresponding to the Wolf algorithm is shown as the 
dashed line. The full line belongs to the new variant. 



Fig. (3) show scatter plots of all three local exponents XI (t), X2 (t), ^3 (t) that have been 
computed for the Rossler attractors. Again, the two parts on the left-hand side show the results for the 
new algorithm and those corresponding to the Wolf algorithm on the right. 

Fig(4) shows the Renyi spectra computed from the point clouds of Figs.( 3.) The dashed lines 
denote for the Wolf algorithm, as before. Fig. (4a) corresponds to Fig.(3) and Fig.( 4b ) the curve 
belonging to the new one, can be observed for small values in three cases (Figs. 2 c, and 4b). In 
general, the calculation of the fractal dimension is less robust (which is between the information 
dimension and the capacity dimension), as discussed in [Kantz andSchreiber,2002]. Systematic errors 
have to be taken into account in those cases. There are some cases where the dashed line 
(corresponding to Wolfs algorithm) increases as a function of q (Figs. 2 b, 2c, and 4a) which indicate 
systematic errors [Ground,2005] . 
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5-Attractor and Bifuraction of Rossler Systems : 

Let us start by briefly describing two typical solutions to the Rossler system [Rossler, 1976] 
readingas (4) where (a, b, c) are the bifurcation parameters. The Rossler system has two fixed points 
given by : 
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x=± ( c ± Vc2 - 4ab)/2 

y = ±(c ± Vc2 - 4ab)/2a (5) 

z=± (c ± Vc2 - 4ab)/2a 



For a = 0.432, b = 2 and c = 4, the Rossler system has a chaotic attractor for solution (Fig. 5a). 
According to Farmer et al. [Farmer and Crutchfield andFroeling and Pachard ,1980], we designate this 
attractor as the spiral attractor. This attractor is characterized by a first-return map to the Poincare 
section. For three-dimensional systems such a section is defined by the plane : 
P = {(yn, zn) € R2|xn = x-, xn > 0} .... (6) 

Thus, the map is constituted by an increasing monotonic branch and a decreasing branch separated 
by the critical point located at the maximum (Fig. 5b). The critical point defines the generating partition 
of the attractor which allows the encoding of all periodic orbits embedded within the attractor [ 
Letellier and Dutertre and Maheu,1995] The increasing branch is close to the bisecting line and, 
consequently, the symbolic dynamics is almost complete. A two-symbol symbolic dynamics [ 
Devaney ,] is complete when all periodic orbits which can be encoded with these two symbols are 
solutions to the Rossler system. Thus, for a=0.432, most of periodic orbits encoded with two symbols 
are embedded within the attractor generated by the Rossler system. 

When the bifurcation parameter a is increased, new periodic orbits are created and the chaotic 
attractor increases in size (Fig. 6b). The corresponding first-return map is constituted by more than two 
branches and, for a = 0.556, up to eleven monotonous branches may been identified [Letellier and 
Dutertre and Maheu,1995]. The corresponding attractor is designated as the funnel attractor [Farmer 
and Crutchfield andFroeling and Pachard ,1980]. For a greater than 0.556, there is metastable chaos, 
that is the trajectory visits the neighborhood of the unstable periodic orbits solution to the Rossler 
attractor before being ejected to infinity [Letellier and Dutertre and Maheu,1995]. The dynamics of 
the Rossler system can therefore be investigated for a < 0.556, b and c remaining constant. 

A bifurcation diagram synthesizes the evolution of the dynamics under the change of the 
bifurcation parameter a (Fig7). The bifurcation parameter a is varied over the interval [0.432, 0.556]. It 
will be shown that quite a similar bifurcation diagram is obtained when the discrimination time step h 
of the discrimination of the Rossler system is increased. 
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Fig. (5). Spiral attractor generated by the Rossler system (5) with the bifurcation parameters (a, 

b y c) = (0.432, 2, 4). 
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Fig.( 6 ). Funnel attractor generated by the Rossler system (5) with the bifurcation parameters (a, 

b, c) = (0.556, 2, 4). 
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Fig. (7). Bifurcation diagram vs. the bifurcation parameter a of the Rossler system (5). Part (a) 

corresponds to values smaller than 0.432 here used 
as a reference and (b) for values larger than this reference. 
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6- Relationship between Sensitivities and initial conditions of Rossler System: 

In this section simulation of chaotic behaviors of Rossler system are depend on definition of 
Gulick . 

The sensitivities are examined by varying the control parameters (a,b,c) . Matlab is used to 
analysis the sensitivities depending on the initial condition . 

(i) (0,0,0.01) and (0.03,0,0.01) with parameter a=10,b=2,c=8 as show in the below figure. 




Fig (8) 

(ii) ( 0,0,0.01) and (0.01,0,0.01) with parameter a=10,b=2,c=8 as show in the below figure 




Fig (9) 



(iii) (0,0,1,0.01) and (0.1,0.1,0.01) with parameter a=0.2,b=0.2,c=6.7 as show in the below figure 




Fig (10) 



(iv) ( 2.05,4,6) and (2, 4,6) with parameter a=0.55,b=2,c=6 as show in the below figure 
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Fig (11) 

(v) ( 0. 1,0,0.3) and (0. 1,0.2,0.3) with parameter a=l,b=2,c=4 as show in the below figure 




Fig (12) 



But, when we change the parameter, we did not get any sensitive to the above initial condition as show 
in the below figures . 





Fig (13) 



Fig (14) 
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